Production of Mesons and Baryons at High Rapidity and High Pt in Proton-Proton 

Collisions at y/s = 200 GeV 
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We present particle spectra for charged hadrons ir , K ± ,p and p from pp collisions at y/s = 200 
GeV measured for the first time at forward rapidities (2.95 and 3.3). The kinematics of these 
measurements are skewed in a way that probes the small momentum fraction in one of the protons 
and large fractions in the other. Large proton to pion ratios are observed at values of transverse 
momentum that extend up to 4 GeV/c, where protons have momenta up to 35 GeV. Next-to- 
leading order perturbative QCD calculations describe the production of pions and kaons well at 
these rapidities, but fail to account for the large proton yields and small p/p ratios. 

PACS numbers: 13.85.Ni, 13.87.Fh, 25.75.Dw, 13.85.Hd, 25.75.-q 



Even though proton-proton collisions are often used as 
a simple hadronic reference system to disentangle nuclear 
effects in p-A or A-A systems, much of that interaction 
remains unknown. Quantum Chromodynamics (QCD), 
the accepted theory of hadronic systems, has had im- 
pressive success in describing many aspects of these col- 
lisions, specially when the underlying partonic processes 
are hard and a perturbative approach is possible. Such 
an approach starts to be applicable at the CERN ISR en- 
ergies ( i/s = 30 — 60 GeV) , but only near mid-rapidity 
[liQ • Measurements of jet and neutral pion production 
from p+p collisions at higher energies at the SppS [3j and 
the TEVATRON [J] show that those systems can be well 
described perturbatively in a wide range of rapidities. 
The RHIC program at BNL has brought the focus back 
to p+p collisions at y/s — 200 GeV with an impressive 
array of probes that range from fully identified hadrons 
in a wide rapidity range, to direct photons, heavy quark 
production and jets. Cross sections measured at RHIC 
have been compared to calculations based on perturba- 
tive QCD (pQCD) where the partonic cross sections are 
calculated up to Next-To-Leading-Order level of accu- 



racy (NLO) [B|, The measured yields of neutral pions 
and direct photons at mid-rapidity @, HI , as well as the 
production of neutral pions at high rapidity are well re- 
produced by the calculations [9|. 

Collider measurements at high rapidities (ybeam = 5.4 
at -y/s = 200 GeV ) offer a window to an interesting re- 
gion of the proton wave function; the kinematics of the 
collision at the partonic level is skewed such that one of 
the protons is probed at small values of fractional longi- 
tudinal momentum x by valence quarks from the other 
proton that have higher values of x. The nucleon inter- 
nal structure has been well measured for a wide range 
of x values, and there are indications that new regimes 
could be reached within the small-x component of the 
wave function where gluons are the most abundant par- 
tons fiol ] . High rapidity measurements in p+p and p+ A 
collisions are thus a useful tool to extend the applica- 
bility of QCD. These data also serve as a reference to 
understand the evolution of the gold wave function with 
rapidity as observed in d+Au collisions, and the forma- 
tion of a dense and opaque colored medium in heavy ion 
collisions [H El . 
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This letter reports the first measurement of unbiased 
differential cross sections for identified charged pions, 
kaons, protons and their anti-particles at high rapidity 
in p+p collisions at s/s = 200 GeV. The comparison of 
these data to a QCD framework furthers the understand- 
ing of the hadron-hadron interaction. 
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FIG. 1: Invariant cross section distributions for pion, kaons 
protons and anti-protons produced in p+p collisions at ^/s = 
200 GeV at rapidities y=2.95 (left panels) and y = 3.3 (right 
panels). In all panels, positive charged particles are shown 
with filled triangles (red in on-line version) and negative ones 
with open circles (blue in on-line version). The errors dis- 
played in these plots are statistical. 

The data were taken with the BRAHMS forward spec- 
trometer set at 4° and 2.3° with respect to the beam 
fl4j |. These settings correspond to y=2.95 and 3.30 re- 
spectively. At these rapidities the geometrical acceptance 
and the particle identification of the experiment allow us 
to reach high values (4-5 GeV/c compared with the 
kinematical limits of 10.4 and 7.4 GeV/c, respectively), 
which is essential to test the ability of NLO pQCD to 
reproduce the measured yields. Measurements at two ra- 
pidities provide a check on rapidly falling cross sections. 

The minimum bias trigger used to normalize these 
measurements was defined with a set of Cherenkov radia- 
tors (CC) placed symmetrically with respect to the nom- 
inal interaction point and covering pseudo-rapidities that 
range in absolute value from 3.26 to 5.25. This trigger re- 
quired that at least one hit is detected in both sides of the 



array. We used events generated with PYTHIA [l~5| and 
a GEANT simulation of the BRAHMS setup to estimate 
that this trigger can access 70 ± 5% of the total inelastic 
proton-proton cross section of 41 mb. The data used for 
this analysis were collected with a spectrometer trigger 
defined by three scintillator hodoscopes, located at 4.6, 
8.8 and 18.9 meters away from the nominal interaction 
point. GEANT simulations of the experiment with input 
events generated by PYTHIA have been used to estimate 
the bias introduced by the CC detectors in the spectrom- 
eter trigger. The deduced correction is equal to 13 ± 7%, 
approximately independent of pt and rapidity, and was 
applied to all spectra. The CC detectors have good tim- 
ing resolution and were also used to define the position 
of the interaction along the beam line with a resolution 
of approximately 2 cm. The present analysis was done 
with charged particles that originated from collisions of 
polarized protons with interaction vertices in the range 
of ±40 cm. The transverse polarization of each beam 
bunch (arranged to point up or down in an alternating 
way) is estimated to have a negligible contribution to the 
inclusive cross-sections presented here: even though the 
bunch intensities for opposite polarizations can differ by 
as much as ~ 10%, the asymmetries of produced particles 
were measured and have small values ( 2-8%) resulting 
in a net effect smaller than 1%. Invariant cross sections 
were extracted in narrow (Ay = 0.1) rapidity bins cen- 
tered at y=2.95 and y = 3.3, respectively. Narrow ra- 
pidity bins are required to reduce the effects of rapidly 
changing cross sections. Each distribution is obtained 
from the merging of up to five magnetic field settings. 
The data are corrected for the spectrometer geometrical 
acceptance, multiple scattering, weak decays and absorp- 
tion in the material along the path of the detected parti- 
cles. Tracking and matching efficiencies for each of the 5 
tracking stations in the spectrometer were calculated by 
constructing full tracks with only 4 stations and evaluat- 
ing the efficiency in the 5*^ station. The overall efficiency 
is about 80-90% and is included in the extraction of the 
cross sections. Particles are identified by our Ring Imag- 
ing Cherenkov detector (RICH) [l6j]. The efficiency of 
this detector has been studied with pions identified with 
a scintillator time-of-flight counter in an overlapping mo- 
mentum range and reaches an upper value of 97%. The 
low momentum part of the proton spectra is measured 
using the RICH in veto mode. 

The momentum resolution of the spectrometers, 5p, 
is calculated from the angular resolution of the track- 
ing detectors. With p in GeV/c, dp/p = 0.0008p at the 
highest field setting. From a comparison of the 3 indepen- 
dent spectrometer measurements and the momentum ex- 
tracted from the ring radius in the RICH we estimate that 
our systematic error on the absolute momentum scale is 
less than 1%. We have identified several sources of pos- 
sible systematic errors and they are listed in Table I . 

Figure 1 shows the invariant cross sections for 
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TABLE I: Estimated systematic errors shown as a per-cent of the measured cross-sections presented in this letter. 
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FIG. 2: Particle ratios versus pr at y=2.95 and 3.3. Top) 
iv~ /tt + , Middle) p/p and Bottom) p/ir + (red circles) and 
p/n~ (blue squares). The shaded rectangles indicate an over- 
all systematic error (17%) estimated for these ratios. The 
dashed line shows an upper limit for the (p + p)/(-K + + tt~) 
ratio from e + e~ collisions. 



7r ± , K ± ,p and p versus pr at y=2.95 and 3.3. The data 
were placed at the average pt value calculated within 
the bin. The upper panels of Fig. [2] show the n~ / n + 
ratio at y—2.95 and 3.3 respectively. Within systematic 
errors both ratios display a falling trend as a function of 
Pt- This may be an indication of the dominance of va- 
lence quark fragmentation at these rapidities and reflects 
the ratio of d to u quarks in the proton. The p/p ratios 
shown in the middle panels are much smaller than unity. 
This is a clear indication that proton and anti-protons 
at these high rapidities are not produced by a common 
mechanism like g — > p fragmentation. The bottom pan- 
els of Fig. [2] show the p/tt + ratio as solid triangles and 
the p/tt~ ratio with solid squares at y=2.95 and y=3.3. 
The dashed line is the maximum value for the measured 
(p + p) / (ir + + 7r~) ratio in e + + e~ annihilation at 91.2 



GeV [17|. The difference between p/ir + andp/7r~ at both 
rapidities is remarkable and caused by an unexpectedly 
large proton yield at these rapidities. Such large yield 
may be related to the mechanism that transfers a con- 



served quantity, the baryon number, from beam to inter- 
mediate rapidities i.e. baryon number transport. What 
remains as an open question is how a mechanism that is 
thought to be mainly restricted to the longitudinal com- 
ponent of the momentum, gives these protons such high 
Pt- 

We have compared our measured differential cross- 
sections with NLO pQCD calculations [IH evaluated 
at equal factorization and renormalization scales, fi 
i i i = fi R = p T . These calculations use the CTEQ6 par- 
ton distribution functions [19j and a modified version of 
the "Knichl-Kramer-Potter" (KKP) set of fragmentation 
functions (FFs) referred to here as mKKP, as well as 
the "Kretzer" (K) set [HI]. The KKP set includes func- 
tions that fragment into the sums tt + +it~ , K + + K~ and 
p + p. Modifications were necessary to obtain functions 
producing the separate charges for both n and K. These 
modifications are based in functional forms supported by 
measurement 22j, and involve the following operations: 
to obtain the FFs producing 7r + , the functions fragment- 
ing favored light quarks u, d into ir° were multiplied by 
(1 + z) ( e.g. Dt = (1 + z)Dl° with Df = \D* + +»") 
where z is the fraction of the parton momentum car- 
ried by the hadron, and D 7 ^ the function fragmenting 
u quarks into positive pions. The functions fragmenting 
unfavored quarks u, d into it are multiplied by 1 — z. The 
same operation is done for n~ , but this time the favored 
quarks are u and d. The FFs of strange quarks and glu- 
ons are left unmodified. Similar modifications were ap- 
plied to obtain FFs into K + and K~ , but this time, the 
starting functions were the ones fragmenting u, u, s, s into 
the sum K + + K~. Figure [3] shows that the agreement 
between the NLO calculations that include the mKKP 
FFs and the measured pion cross section is remarkable 
(within 20% above 1.5 GeV/c) even though the modi- 
fications applied can be considered as "crude" . Similar 
ood agreement was obtained for neutral pions at y=0 
and at y=3.8 [9] at RHIC. The agreement between 
the calculated and the measured kaon cross-sections is 
equally good. The difference between the mKKP and K 
parametrizations is driven by higher contributions from 
gluons fragmenting into pions. This difference has been 
identified as an indication that the gg and gq processes 
dominate the interactions at mid-rapidity [7|; these re- 
sults indicate that such continues to be the case at high 
rapidity. The calculation that uses the K set underesti- 
mates the pion yields by a factor of ~ 2 at all values of 
Pt while for positive kaons, the agreement is good at low 
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FIG. 3: Top) Comparison of invariant cross sections for tt~ , K + , p and p at y=2.95 and NLO calculations with factorization 
and renormalization scales set equal to pr- The mKKP set of fragmentation functions (solid red line on-line) produce the best 
agreement with the ir~ and K + data. The p and p are compared with the calculation using the AKK set divided by 2 (dashed 
red line in on-line version), see text for details. Bottom) Relative differences between data and calculations. The top smooth 
curves show the effect of setting /j, — 2pr and the bottom curves p = l/2px- For the baryons the (red) filled triangles show p 
data vs the AKK/2 set. 



momentum but deteriorates at higher momenta. 

An updated version of FFs that we refer to as the 
"Albino, Kniehl and Kramer" (AKK) set has been ex- 
tracted from more data made available recently [23(. It 
reproduces well the p + p distributions measured at mid- 
rapidity by the STAR collaboration 24j. At high rapid- 
ity, the contribution from gluons fragmenting into p or p 
is dominant in this new set of FFs (> 80% for px < 5 
GeV/c [IH), and the calculated cross sections for both 
particles consequently have nearly the same magnitude. 
We thus compare the measured cross sections for p and 
p to the NLO calculation using the AKK FFs divided by 
2 in the right-most panel of Fig. [3] The calculation is 
close to the measured p cross section but it is almost an 
order of magnitude higher than the measured p distribu- 
tion. We conclude that the AKK FFs cannot be used to 
describe baryon yields at high rapidity because they fail 
to reproduce the measured abundance of p with respect 
to p. We have ruled out the use of the standard KKP 
FFs because they produce p + p cross sections that are 
smaller by a factor of ~ 10 compared to the measurement 
(not shown). 

In summary, unbiased invariant cross sections of identi- 



fied charged particles as function of pr were measured at 
high rapidity in p+p collisions at y/s = 200 GeV. NLO 
pQCD calculations reproduce reasonably well the pro- 
duced particle (pions and kaons) distributions but p and 
p cannot simultaneously be described well by any of the 
available FFs. These results may show a limitation of 
the factorized description of p+p cross sections, perhaps 
because it does not include the effects of baryon number 
transport that, as the data suggest, may extend to high 
Pt- These measurements bring additional insight into 
the hadron-hadron interaction and its description in the 
context of QCD; they are as well instrumental in con- 
straining phenomenological descriptions of that system. 
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